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ABSTRACT 

Quasars represent a brief phase in the life-cycle of most massive galaxies, but the evolutionary connection be- 
tween central black holes and their host galaxies remains unclear. While quasars are active and shining brighter 
than the Compton-limit luminosity, their radiation heats the surrounding medium to the Compton temperature, 
forming Compton spheres extending to the Stromgren radius of Fe 26+ /He 2+ . After the quasars shut off, their 
"afterglow" can be detected through three signatures: (1) an extended X-ray envelope, with a characteristic tem- 
perature of ~ 3 x 10 7 K; (2) Lya and Ly/3 lines and the K-edge of Fe 26+ ; and (3) nebulosity from hydrogen and 
helium recombination emission lines. We discuss the possibility of detecting these signatures using Chandra, the 
planned XEUS mission, and ground-based optical telescopes. The luminosity and size of quasar afterglows can be 
used to constrain the lifetime of quasars. 

Subject headings: black hole physics — cosmology: theory — galaxies: formation — quasar: general 
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1. INTRODUCTION 

The neutral hydrogen bounded by dark matter halos at high 
redshift has detectable extended Lya "fuzz" powered by pho- 
toionization from central quasars (Haiman & Rees 2001). The 
detection of this signal may impose strong constraints on quasar 
formation during galaxy assembly. The lifetime of quasars re- 
mains uncertain, but they are generally thought to be short- 
lived, with estimated ages of ~ 10 7 to 10 9 years (Jakobsen et 
al. 2003; Martini & Schneider 2003). Detecting the infalling 
gas, heated by the quasar, as Lya fuzz can help to constrain the 
quasar's lifetime. If undetected, it may imply that quasar activ- 
ity only appears when the gas has settled into a thin disk or has 
already been consumed mostly by star formation. The prelimi- 
nary evidence for extended Lya emission in radio-quiet quasars 
already suggests that quasar activity influences galaxy assem- 
bly (Bergeron et al. 1999; Fardal et al. 2001; Bunker et al. 
2003). The spectrum of the quasar Q1205-30 unambiguously 
shows an extended hydrogen Lya glow formed by gas falling 
into the ionizing cone of the quasar (Weidinger et al. 2004). 

What is the physical appearance of a quasar turning o/f within 
an assembling galaxy? This Letter investigates the afterglow 
properties of ultraluminous, radio-quiet quasars. 

2. COMPTON SPHERE AND COMPTON LIMIT 

While the quasar is active, it provides a source of energy to 
heat infalling material (Binney & Tabor 1995; Ciotti & Ostriker 
2001; Brighenti & Mathews 2002). Hot gas in elliptical galax- 
ies has been extensively studied for three decades since Ein- 
stein, but the physics of its origin and evolution still are not 
completely understood (see the review by Mathews & Brighenti 
2003). It has been realized that gas infalling into the dark mat- 
ter halo of a nascent galaxy (Fall & Rees 1985) undergoes ther- 
mal instability and forms a two-phase medium that would ex- 
hibit very faint emission features (Haiman, Spaans & Quataert 
2000). In this stage, pressure balance holds between the cold 
(T c ss 10 4 K) and hot (7], ss 10 6 K) medium. The clumping fac- 
tor of the two-phase medium is C = (n 2 ,)/ («t>) 2 = (Th/T c ) ~ 



10 4 , where «b is the baryon number density. The irradiation 
by the quasar light greatly changes the two-phase structure of 
the infalling gas. We assume the abundance of hydrogen, he- 
lium and iron to be An = 0.7, _4 H e = 0.3 and _4 Fe = 5.0 x 10" 5 , 
respectively. The thermal state of the photoionized gas is deter- 
mined by the ionization parameter, defined by 

cnkT ' ^ ^ 

where n and T are the number density and temperature of the 
gas and F is the flux of the ionizing source (Rrolik, McKee & 
Tarter 1981). Compton heating and cooling balance at a tem- 
perature of 7c = e/Ak, where e is the mean energy of the ioniz- 
ing source and k is Boltzmann's constant. For a typical quasar 
spectrum (Elvis et al. 1994), the mean photon energy is e ss 10 
keV. The gas will be heated to Tq ~ 3.0 x 10 7 K (see the de- 
tailed calculation of Sazonov et al. 2004), which is one order 
of magnitude hotter than the typical virial temperature of 0.15 
keV for the hot gas (Fall & Rees 1985). When the ionization 
parameter of the ionized gas reaches 

I --3/2 



3 > 3 C = 6.1?! 



(2) 



where l\ = e/lOkeV, the gas will be fully heated to the Comp- 
ton temperature (Rrolik, McKee & Tarter 1981). The timescale 
of Compton heating is given by f c eat = 3m e c 2 kT c /Sa T eF, where 
m e is the electron mass and aj is the Thomson scattering cross 
section. This timescale is derived by dividing the number of 
scatterings needed to reach the Compton temperature by the 
scattering rate. We assume a quasar lifetime of ?q = TQfs a i P , 
with ?sai P = 4.6 x 10 8 yr being the Salpeter timescale (Salpeter 
1964). Setting f c eat = f Q , with L 46 = L/10 46 erg s _1 , we have the 
Compton radius 



R c = 3.2L% 2 T ^ 2 kpc 



(3) 

which, along with the condition 5 = 3c, gives the Compton 
density 

cm -3 . (4) 



«c : 



Lixio-V^ep 
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FIG. 1. — Illustration of the afterglow from an ultraluminous quasar. The 
Compton radius, Rq, extends to the Stromgren radius, R$, during the quasar's 
lifetime when its luminosity exceeds the Compton limit luminosity. The after- 
glow arises from the Compton sphere after the quasar is extinct. 

The Compton radius depicts the size of the region in which 
the gas is heated by the quasar to the Compton temperature dur- 
ing the timescale fQ. The Compton density is the maximum 
density allowed within Rc by Compton heating during the life- 
time of the quasar. It is interesting to note that nc is independent 
of the quasar luminosity. Beyond Rc, the gas has an ionization 
parameter of 5 < Sc and hence is in a two-phase state. For a 
medium with a density lower than tic, the size of the region with 
S > Sc extends beyond Rc- However, in such a case Comp- 
ton heating is sufficiently inefficient that the gas is in a time- 
dependent two-phase state. We do not consider this regime for 
ultraluminous quasars. 

The radius of the Stromgren sphere of Fe 26+ ions is given by 



R 



FeXXVI , 



4-Ka Fe A Fe C(n 2 b ) 



1/3 



- 5 7 L 1/3 C" 1/3 



(4) 



-4 



kpc, 
(5) 



26* 



where C 4 = C/10 4 , a Fe = 9.06 x 10" 12 cm 3 s" 1 is the Fe 
recombination coefficient in Seaton (1959) evaluated at T e = 
3 x 10 7 K, and (n|}-4 = (n|}/10~ 4 cm" 6 is the volume-averaged 
mean square density of the cold gas. The rate of ioniz- 
ing photons, Nioa, is estimated from the typical spectrum of 
quasars (Elvis et al. 1994). The radius of the He 2+ Stromgren 
sphere Rf m = 10.3 L%% 1/3 (nf)^ 3 kpc, and R^/Rf 1 » 
0.3; we assume ionization energies of Enem = 54.4 eV and 
Emi = 13.6 eV and an ionizing spectrum /„ oc t/~ 1,3 . The ra- 
dius of the ionized region is given by (Cen & Haiman 2000) 

Hem_^M t _lA„A__ /„2\l/2^/3 = r l/3 L l/3 /ii2 \l/6 



■46 ( n l)-4 k P C < 



and ^n IU /^Sn = 1-2- As the ionizing luminosity increases, the 
Compton radius extends to the radius of the helium Stromgren 
sphere. The entire Stromgren sphere will be fully heated to the 
Compton temperature when R^ eXXV1 = R c , at which point the 
iron Compton-limit luminosity 

(6) 



L p c e = 3.0 x 10 4 ' T-'C- A l {niY_i erg s 
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and the iron Compton-limit radius 



R% = 11.7^(4) 



2 U kpc. 



(7) 



FIG. 2. — The X-ray spectrum of the Fe 26+ Compton sphere. The green line 
represents free-free emission, the red line free-bound (radiative recombina- 
tion) emission, and the black line the total emission from the Compton sphere. 
Iron Lya A 1.78 and Ly/3 A 1.51 lines, with emissivities (Raymond & Smith 
1977) for ~ 3 X 10 7 K and solar iron abundances, are labeled (AE/E = 0.02 
for XEUS). The Lyo and Ly/3 lines will have equivalent widths of 12.0 keV 
and 6.6 keV for solar iron abundances, and should be detectable with XEUS. 

When L < L c e , the Compton sphere is smaller than the 
Stromgren sphere. Since the medium in the region between the 
Stromgren radius and the Compton radius cannot fully reach 
the Compton temperature, the features of the afterglow arising 
from it could be complicated. We are currently interested in the 
case of L > L c e , namely, the ultraluminous quasars. 

The energy stored within the Compton sphere, ^nR c n e kT « 

4 x 10 58 n_2 erg, where T = 3.0 x 10 7 K, «_ 2 = n e /10" 2 cm -3 , and 
Rc = 20 kpc, is only a small fraction (0.02%) of the total energy 
released (2 x 10 62 erg) during the activity of a quasar with a 
black hole mass of 1O 9 M0 and an accretion efficiency of 0.1. 
The recombination timescale is {a Fe n e )"' w 0.3 x 10 7 «l2 Y 1 "- F° r 
quasars surrounded by a very dilute medium, this timescale is 
too long for the ions to recombine within the Hubble time at 
the relevant epoch. The Stromgren radius is then never reached, 
and the hot gas would emit no recombination lines. This situa- 
tion might be relevant for some isolated, exceptionally gas-poor 
elliptical galaxies. If the recombination timescale is longer than 
the lifetime of quasars, the discussion on the Stromgren radius 
is also inapplicable. 



3. OBSERVATIONAL SIGNATURES 

An afterglow will be seen from the Compton sphere after 
the quasar is extinct. Here we estimate the expected signature 
of the quasar afterglow from the hot gas inside the Compton 
sphere that is in hydrostatic equilibrium with the gravitational 
potential of the dark matter halo (Fig. 1). The isothermal den- 
sity profile of the hot gas is given by (Makino et al. 1998) 



p g = poe-P (1+R/R, 



f Rs/R , where p = 1.96 x 10" 26 gem" 3 , 

fl/3 



= 3.Q3M n T^\ the characteristic radius R s = 44.25M^ J kpc, 
T 7 = T/10 7 K, and M n = M halo /10 13 M . We use a concentra- 
tion parameter of cq = 5 and A c = 18tt 2 for the dark matter halo. 
The gas mass fraction is only ~ 0.005(r es /10~ 3 ) of that of the 
dark matter halo, where r es is the Thomson scattering depth of 
the hot gas. We note that the density of the hot gas is lower than 
the Compton density. 

The cooling timescale of the fully ionized gas inside the 
Compton sphere is governed mainly by free-free radiation. 
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FIG. 3. — Expected fluxes of hydrogen (Lya A1216) and helium (He II 
AA1640, 4686) recombination lines. The solid and dotted lines show the 
predictions for an electron temperature of 10 4 K and 2 X 10 4 K, respec- 
tively. We assume a ACDM cosmology with f2 m = 0.3, Qjy = 0.7, and 
H = 7510ns- 1 Mpc" 1 . 

The cooling timescale at the center of the halo is fff w 

1.0 x l0 9 Tj^ 2 nz\ yr, and it will be considerably longer in the 
outskirts. This leads to the appearance of a long-lived, ex- 
tended X-ray-emitting envelope, characterized by a luminosity 
Lff w 1.3 x 10 42 erg s" 1 for T = 10 7 K, Rc = 2Q kpc and the above 
density profile, peaking at an energy of a few keV. Since the de- 
gree of ionization is very high during the active phase of the 
quasar, we assume that iron is fully ionized. Figure 2 shows the 
predicted X-ray signature from such a model, where we have 
computed the integrated emission from the Compton sphere of 
a progenitor quasar shining at a fiducial Compton luminosity 
of 10 42 erg s" 1 , including contributions from the K and L edges 
of Fe 25+ free-bound emission from hydrogen-like ions (Seaton 
1959) as well as the Lya and Ly/3 lines from Fe 26+ . For redshifts 
Z < 0.2, the expected flux of > 1.1 x 10~ 14 erg s" 1 cm" 2 will 
be detectable by Chandra/ AGS with an exposure time of 100 
ks, and the angular radius of the Compton sphere, 9 w 29.5", 
should be resolvable. However, the effective area of ACIS 
(~ 80 cm 2 above a few keV) is too small to be able to de- 
tect the Lya and Ly/3 lines of Fe 26+ . The future XEUS mis- 
sion, with an expected angular resolution of 2" at a flux level 
of 10~ 17 erg s" 1 cm" 2 in the energy range 0.05-30 keV, will 
be able to detect the Compton sphere and Lya and Ly/3 emis- 
sion at redshifts lower than 1. Higher-redshift afterglows are 
too faint to be detected by Chandra, whereas XEUS will have 
insufficient resolution. 

Lastly, we consider the signal from the region between the 
Compton radius and the ionization radius (Fig. 1), which arises 
from hydrogen and helium photoionized by the quasar. Af- 
ter the quasar quenches, the photoionized gas will radiate re- 
combination lines of hydrogen and helium. We calculate the 
fluxes of hydrogen Lya and He II AA1640, 4686, using the re- 
combination coefficients in Seaton (1978), by integrating the 
shell between the Compton and ionization radius (Fig. 3). The 
broadening width of the recombination lines due to thermal mo- 
tions, AA/A w (kTc/ZmpC 1 ) 1 ' 1 w (0.2-3.0) x 10" 3 , where Z 
is the charge of the ions. This width can be readily resolved 
by ground-based optical telescopes, the largest of which should 
also have the requisite sensitivity and spatial resolution to im- 
age the line-emitting region. We note that the emission-line 
regions considered here are different from the extended Lya 



nebulae discovered in some star-forming galaxies (Steidel et al. 
2000). The latter are powered by OB stars, which are not ener- 
getic enough to ionize He. 

If quasar afterglows can be detected, either the luminosity 
or the size of the Compton sphere (Equations 6 and 7) can 
be used to estimate tq, the lifetime of the progenitor quasar. 
In practice, the luminosity would be much better constrained 
than the size, as the latter would be only marginally resolved 
even with Chandra and depends critically on the detailed sur- 
face brightness distribution. While the hot gas normally associ- 
ated with individual elliptical galaxies (O' Sullivan et al. 2003) 
rarely has temperatures exceeding 1 keV, galaxies situated in 
groups or denser environments would be surrounded by hotter 
gas, which at large distances would be difficult to distinguish 
from any component intrinsic to an afterglow. To minimize 
this source of confusion, the local environment of any poten- 
tial targets needs to be carefully scrutinized. With typical X-ray 
temperatures of a few keV, X-ray binaries present another pos- 
sible source of contamination, but their integrated X-ray emis- 
sion (Fabbiano 1989), even in large galaxies, is significantly 
less than 10 42 erg s" 1 . 

If quasars shine near the Eddington limit, the Compton- 
limit luminosity corresponds to a black hole mass Mbh ~ 
2 x 10 9 Mq. If such massive black holes are common among 
quasars (Netzer 2003; Vestergaard 2004), their afterglows 
should not be difficult to detect. From the Mbh - o relation 
(Tremaine et al. 2002), we expect the host galaxies to have 
stellar velocity dispersions of a w 400 km s" 1 . Thus, future 
searches of quasar afterglows should target the most massive 
giant ellipticals. Such galaxies are rare, but not exceedingly 
so. From the local velocity dispersion function of early-type 
galaxies (Sheth et al. 2003), the space density of galaxies with 
a « 400 km s" 1 is - 10" 6 Mpc" 3 . 

For the more abundant but less powerful quasars, with lu- 
minosities below the Compton-limit luminosity, the Compton 
radius cannot match the Stromgren radius and it will be signif- 
icantly smaller than the Compton-limit radius. The medium in 
the transition zone between the Compton and Stromgren radii is 
in a time-dependent two-phase state, since it is not able to reach 
an equilibrium state within the lifetime of the quasar. We specu- 
late that these lower-luminosity sources would have optical and 
ultraviolet emission features quite similar to those of ultralumi- 
nous quasars, but that the X-ray spectrum should be dominated 
by the Lya and Ly/3 lines of lighter ions, such as Ne X or O VIII. 
Detailed calculations for this case are in progress. 

4. SUMMARY 

In conclusion, we show that the afterglows of ultraluminous 
quasars form relic Compton spheres that could be detected by 
Chandra, XEUS, and ground-based optical telescopes. The lu- 
minosities and sizes of the afterglows can be used to estimate 
the lifetime of quasars, a key constraint on quasar and galaxy 
evolution. The extended X-ray envelope recently detected 
around the isolated elliptical galaxy NGC 4555 (O' Sullivan & 
Ponman 2004) could be an example of the afterglows discussed 
in the present paper. It would be worthwhile to search in this 
object for the predicted extended emission from hydrogen and 
helium recombination lines. 
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